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Abstract 



I compare the lattice calculation of Fukugita et al. for the pion S'-wave scattering 
lengths to the predictions of Chiral Perturbation Theory to two loop accuracy. I 
find good agreement, despite the use of the quenched approximation in the lattice 
calculation. 



PACS numbers: 12.39.Fe, 12.38.Gc 



In this letter I will present the comparison of two results which are already 
available in the literature: the calculation of the pion S'-wave scattering lengths in 
quenched lattice QCD |]T| (see also [^) and in Chiral Perturbation Theory (CHPT) 
to two loops 0. Admittedly, the two results need not be the same: the physical 
content of the two calculations is different since QCD contains sea quark contri- 
butions, whereas the quenched lattice calculation does not. Let me disregard this 
complication for the moment, I will come to it later. 

There is an important reason for attempting a comparison, which however con- 
cerns more future developments than present day calculations. As it is well known 
CHPT provides an efficient technique for implementing the constraints of the chiral 
symmetry of QCD on the Green functions. Remarkably, at low energy this is enough 
to make predictions. This approach uses only the symmetry properties of QCD: the 
dynamics of quarks and gluons is totally ignored and shows up only in the values 
of some unknown low energy constants, which we usually learn from experiments. 
The major challenge for bridging the gap between QCD and its low energy effective 
theory is to reproduce the value of these low energy constants starting from the 
QCD Lagrangian. The best candidate for doing this is lattice QCD, although very 
little has been done up to know - the only attempt in this direction of which I am 
aware is the work by Myint and Rebbi |Q. 

The comparison of lattice QCD and CHPT for quantities like the pion S'-wave 
scattering lengths is very interesting from this point of view: finding agreement 
means that we understand the value of these low energy constants - as extracted 
from experiments - on the basis of the QCD Lagrangian. Of course the full accom- 
plishment of this program requires overcoming the quenched approximation. 

The method to calculate scattering lengths on the lattice is due to Liischer 0. 
In principle it allows to extract these quantities in a theoretically very clean manner. 
The method consists in measuring the energy of the fundamental state of two static 
particles inside a spatial box of length L. Liischer has calculated the expansion of 
this energy in inverse powers of the length L, showing that the first three coefficients 
of this expansion (which starts with a term of order L~^) contain powers of the S'- 
wave scattering length. From the measurement of the energy shift due to the finite 
size of the box one can extract the scattering length. The method is very powerful, 
in particular because it allows to directly relate a calculation done in Euclidean 
space to a quantity which is defined in Minkowski space. Moreover it transforms 
what is usually considered a systematic error into an interesting effect from which 
to extract useful information. 

Liischer's formula can be applied only if M^L ^ 1. This condition was satisfied in 
the calculation of Ref. |l|] by using an unphysical, moderately large pion mass. This 
is not a problem for comparing to the CHPT calculation, since in the chiral expansion 
the dependence on light quark masses is fully explicit: given m = 1/2 (m^j + m^), or 
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a value for and F^r, one can calculate the corresponding value for the scattering 
lengths. The only requirement for the chiral expansion to be justified is that the pion 
mass be small with respect to the typical mass scale of QCD, Mp say. This condition 
is satisfied in the case of staggered fermions, where Fukugita et al. had M^/Mp ~ 0.3. 
On the contrary, we cannot compare to the Wilson fermions calculations, where 
M^/Mp ~ 0.7. 

Let us discuss here in detail how the CHPT expressions have to be calculated in 
the case of unphysical quark masses. For the purpose of this discussion I will only 
consider the one loop case. The extension to two loops is straightforward. At one 
one loop Oq is given by [0: 




y(/U/i) + /^(/^)) + y/S(/i) + 2/K/^) 
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This expression is scale independent: the scale dependence of the low energy con- 
stants /[(/i) compensates the one due to the chiral log. 

To compare to the lattice result I use the ratio M^/F^ as calculated on the 
lattice. The ^[(/i) do not depend on the light quark masses and I use for them the 
values given in Ref. [Q. The last thing to be evaluated is the chiral log. This requires 
to give the value of M^^ in physical units. To make the connection between lattice 
and physical units I use the p mass. This depends on the light quark masses and can 
be expressed as a power series. Using an ansatz known as the "quark counting rule" 
(see e.g. Ref. |^), I consider only the first two terms in the expansion, according to: 

Mp = My + 2m + 0(m^/2) , (2) 

where My stands for the p mass in the chiral limit, for which I use the value My = 
760 MeV. According to this ansatz, for low values of M^^/Mp the Mp dependence 
on the light quark masses makes a rather small effect, and the pion mass in physical 
units, at MjMp = 0.33 turns out to be M^ = 260 MeV (using also B = 2My [§, 
where B = — (0|gg|0)/-F^). It should be noted that the higher the value of the pion 
mass, the higher is the relative importance of the low energy constants with respect 
to the chiral log (at fixed p). This means that on the lattice one can in principle 
increase the sensitivity of a given quantity to the low energy constants, and improve 
the accuracy of their determination. 

The same procedure has to be repeated in the two loop case. There we have six 
new constants appearing, the r[(/i), z = 1, . . . 6 whose value has been estimated via 
resonance saturation in Ref. [^. Here I use the same estimate. The effect of these 
new constants is small. 
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Figure 1: Comparison of the / = pion scattering length as calculated on the 
lattice by Fukugita et al. DH and the CHPT expansion. The data point with the 
square is the experimental measurement The two lattice points are obtained 
with staggered fermions. The circle refers to the calculation done without gauge 
fixing, and the triangle to the Coulomb gauge. 

The curves for the dimensionless quantities 32F^aQ/M^ as functions of M^^/Mp 
are shown in Figs. 1, 2 for J = 0, 2 respectively [|ll[|. The straight solid line corre- 



sponds to the Current Algebra {i.e. CHPT to tree level) prediction [12|, whereas 
the dotted curve to the one loop CHPT calculation. The three dashed curves are 
all calculated with the two loop expression for the scattering length. The central 
one corresponds to the values of the constants used in Ref. @], whereas the other 
two are obtained by varying the low energy constant Zgl^p) by ±2 10~^: this is a 
generous error for this particular constant, and is meant to give a rough estimate of 
the overall uncertainty as coming also from the other low energy constants. 



In Figs. 1 and 2 I have shown also the available experimental data points |T0 
and the lattice points calculated by Fukugita et al. [Q. In the J = case the two 
loop CHPT curve is somewhat higher than the lattice points. In addition, if one 
looks at the relative contribution of the two loops with respect to the one loop and 
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Figure 2: Comparison of the 1 = 2 pion scattering length as calculated on the lattice 
by Fukugita et al. [||] and the CHPT expansion. The legend for the data points is 
the same as in Fig. 1; the experimental datum is from Ref. 

tree level, one can see that the series is not converging very rapidly, and one may 
expect positive sizeable contributions from higher orders. In the 1 = 2 case the one 
and two loop corrections to the Current Algebra value are rather small, and the 
lattice calculation is in agreement with all of them. All in all, the lattice calculation 
and CHPT to two loops show deviations from Current Algebra which agree in sign 
and size surprisingly well, for both isospin cases. This seems to indicate that the 
effects due to quenching are of the same size, or smaller than the statistical errors 
of the lattice calculations. 

Finally, I come back to the question of a priori estimates of the effect of quench- 
ing. In a recent article Bernard and Golterman [|T3|, have analyzed the modifications 



of Liischer's formula in the presence of quenching. For this purpose they used the 
method called quenched CHPT In this approach one is able to obtain two 



types of results: first, one can see whether in a specific quantity the quenched ap- 
proximation removes the chiral logs (in Ref. |T^ they showed that in M^^ and 



to one loop there are none); and secondly one can estimate the effects due to the 
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?7-singlet propagator which in this approximation is ill-defined (because it contains 
a double pole). These results show that the quenched approximation introduces a 
qualitative change in approaching the chiral limit. For finite quark masses, how- 
ever, it is difficult to estimate reliably the size of the change: in CHPT the chiral 
logs come together with the low energy constants that remove the scale dependence 
(see Eq. (|l])), and the splitting between the two is arbitrary. At most one can try 
to estimate the part of the effect which is due to the double pole in the ?7-singlet 



propagator (see Ref. [T^ for details). This is what Bernard and Golterman have 



done in the speciffc case of Liischer's formula ||13[. Their result is that the quenched 
approximation modiffes considerably the expansion of the energy shift in inverse 
powers of the box length: in the 1 = case, for example, it introduces terms of 
order l/L"" with n = 0,2. Nevertheless their numerical estimate of these spurious 
effects is not too discouraging: for the staggered fermions calculation of Fukugita et 
al. they estimate a 1% effect in the / = case (because of a strong cancellation), 
and a 20% correction in the 1 = 2 case. 

This a priori estimate suggests that the quenched approximation does not spoil 
completely the results of this lattice calculation, as the comparison proposed here 
also indicates. On the other hand, the use of quenching certainly diminishes the 
importance of the agreement with two loop CHPT. In fact, had this result been 
obtained in full QCD, one could conclude that lattice QCD is able to explain the 
value of the two combinations of low energy constants appearing in the two scattering 
lengths: this would be a very remarkable achievement of lattice QCD. Unfortunately, 
in the present case I cannot go that far, and I have to conclude with the hope that 
there will be soon various improvements on this lattice calculation. For example, 
a repetition of this calculation with a different lattice size would allow one to see 
whether the terms of order l/L" with n = 0,2 show up, and eventually to remove 
them. Moreover it would allow a more reliable extraction of the coefficient of the 
term, by explicitly checking the volume dependence. Hopefully, this and other 
future improvements will tell us whether the agreement observed here contains some 
real physics, or is only accidental. 



Acknowledgements I thank for interesting discussions and/or support J. 
Gasser, P. Hasenfratz, A. Papa, R. Petronzio and R. Urech , and J. Gasser for very 
useful comments on the original manuscript. This work has been supported in part 
by Schweizerischer Nationalfonds and HCM-EEC-Contract No. CHRX-CT920026 
(EURODA$NE). 



References 



6 



[1] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino, and A. Ukawa, Phys. Rev. 
Lett. 71 2387 (1993); Phys. Rev. D 52 3003 (1995). 

[2] M. Guagnelli, E. Marinari, and G. Parisi, Phys. Lett. B240 188 (1990); 
S. R. Sharpe, R. Gupta, and G.W. Kilcup, Nucl. Phys. B383 309 (1992); 
R. Gupta, A. Patel, and S. R. Sharpe, Phys. Rev. D 48 388 (1993). 

[3] J. Gasser, and H. Leutwyler, Ann. of Phys. (N.Y.) 158 142 (1984). 

[4] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, and M. Sainio, Phys. Lett. B 
374 210 (1996). 

[5] S. Myint, and C. Rebbi, NucLPhys. B421 241 (1994). 

[6] M. Liischer, Comm. Math. Phys. 105 153 (1986). 

[7] J. Gasser, and H. Leutwyler, Phys. Lett. 125B 325 (1983). 

[8] J. Gasser, and H. Leutwyler, Phys. Rep. C87 77 (1982). 

[9] L. Rosselet et al., Phys. Rev. D 15 574 (1977). 

[10] M.M. Nagels et al., Nucl. Phys. B147 189 (1979). 

[11] To obtain the curves I had to vary the ratio M^^jF.^ continuously, letting it 
go through 2.24 at M^jMp = 0.18 (physical value), and 4.8 at M^/Mp = 0.33 
(lattice value, see Ref. ll|). How I interpolated between these values is irrelevant 
for the present discussion. 

[12] S. Weinberg, Phys. Rev. Lett. 17 616 (1966). 

[13] G.W. Bernard, and M.F.L. Golterman, Phys. Rev. D 53 476 (1996). 
[14] G.W. Bernard, and M.F.L. Golterman, Phys. Rev. D 46 853 (1992). 



7 



